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ABSTRACT 



> 

m 

, The Local Group galaxies constitute a fundamental step in the definition of cosmic distance scale. Therefore, obtaining accurate distance 
«^ ■ determinations to the galaxies in the Local Group, and notably to the Andromeda Galaxy (M31), is essential to determining the age and 
evolution of the Universe. With this ultimate goal in mind, we started a project to use eclipsing binaries as distance indicators to M31. Eclipsing 
, binaries have been proved to yield direct and precise distances that are essentially assumption free. To do so, high-quality photometric and 
spectroscopic data are needed. As a first step in the project, broad band photometry (in Johnson B and V) has been obtained in a region (34' x34') 
at the North-Eastern quadrant of the galaxy over 5 years. The data, containing more than 250 observations per filter, have been reduced by 
means of the so-called difference image analysis technique and the DAOPHOT program. A catalog with 236 238 objects with photometry in 
Q ' both B and V passbands has been obtained. The catalog is the deepest (V < 25.5 mag) obtained so far in the studied region and contains 3 964 
identified variable stars, with 437 eclipsing binaries and 416 Cepheids. The most suitable eclipsing binary candidates for distance determination 
CZ2 ■ have been selected according to their brightness and from the modelling of the obtained light curves. The resulting sample includes 24 targets 



^ ■ 

Or 



with photometric errors around 0.01 mag. Detailed analysis (including spectroscopy) of some 5-10 of these eclipsing systems should result in 
a distance determination to M31 with a relative uncertainty of 2-3% and essentially free from systematic errors, thus representing the most 
accurate and reliable determination to date. 
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1. Introduction 

Extragalactic distance determinations depend largely on the 
calibration of several distance indicators (e.g., Cepheids, su- 
pernovae, star clusters, etc) on galaxies of the Local Group with 
known distances. The Large Magellanic Cloud (LMC) has tra- 
ditionally been used as a first rung for extragalactic distance 
determinations. However, even taking into account that some 
recent results seem to converg e to an LMC distance modulus 
of (m - M) = 18.50 + 0.02 jAlvesll2004 . its low metallic- 
ity and a possible line-of-sight extension have posed serious 
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difficulties to the accurate calibration of several distance in- 
dicators. The Andromeda galaxy (M31), on the contrary, has 
a more simple geometry and a metallicity more similar to the 
Mil ky Way and other galax ies used for distance estimation (see, 
e.g.. iFreedman et aflll994l) . Therefore, although stars in M31 
are about six magnitudes fainter than those in LMC, the partic- 
ular characteristics of this galaxy make it a promising first step 
of the cosmic distance scale dClementini et al.l200lh . 

Given the importance of M3 1 as an anchor for the extra- 
galactic distance scale, many studies have provided distance 
determinations to M31 using a wide range of methods. A com- 
prehensive list of distance determinations to M31, with explicit 
errors, are shown in Tabled As can be seen, the values listed 
are in the range (m - M)$ = 24.0 - 24.6 mag. Although the 
weighted standard deviation is of -4%, most of the distance de- 
terminations in Table^rely on previous calibrations using stars 
in the Milky Way or the Magellanic Clouds. As a consequence 
of this, a large number of subsequent distance determinations, 
based only on new recalibrations, can be found in the litera- 
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ture. These are not included in Tabled Therefore, a direct and 
precise distance determination to M3 1 is of central importance 
since this would permit the use of all the stellar populations in 
the galaxy as standard candles. 

Table 1. Distance determinations to M31 as presented in the 
references. Values resulting of posterior calibrations and dis- 
tance moduli without extinction corrections are not included. 



Method 


(m - M) 


Distance 


Reference 




[mag] 


[kpc] 




Cepheids 


24.20+0.14 


690+40 


[1] 


1 lp ot the KOB 


Z4.4U±0.ZJ 


/ou+yu 


|2J 


Cepheids 


24.26±0.08 


710±30 


[3] 


RR Lyrae 


24.34±0.15 


740±50 


[4] 


Novae 


24.27±0.20 


710±70 


[5] 


Cepheids 


24.33+0.12 


730±40 


[6] 


Cepheids 


24.41 ±0.09 


760±30 


[6] 


Cepheids 


24.58±0.12 


820±50 


[6] 


Carbon-rich stars 


24.45±0.15 


780±50 


[7] 


Cepheids 


24.38±0.05 


752+17 


[8] 


Carbon-rich stars 


24.36±0.03 


745±10 


[8] 


Glob. Clus. Lum. Func. 


24.03+0.23 


640±70 


[9] 


Red Giant Branch 


24.47±0.07 


780±30 


[10] 


Red Clump 


24.47±0.06 


780+20 


[11] 


Red Giant Branch 


24.47±0.12 


780+40 


[12] 


Cepheids 


24.49±0.11 


790±40 


[13] 


RR Lyrae 


24.50+0.11 


790±40 


[14] 


Tip of the RGB 


24.47±0.07 


785+25 


[15] 


Eclipsing binary 


24.44+0.12 


770+40 


[16] 


Mean & std. deviation 


24.39+0.08 


750+30 





[1J 

[3] 



Baade Xs^^fT^l : r21:lMould & KristianlfT 
Welch etalJll98d): Pritche t & van den Ben 
[5]: Capaccioli et alJ <19 89>: [6]: F reedman & Madori 
[7] 
[9] 



Richer et alJll990|): r81: | Brewer et alJ<199f 
Ostriker & GnedinN 19971: [1 01:lHollandHl99?^ 




[11]: Stanek & GarnavichNl998t): ri21:lDurrell et alJfcOOll) : 



[13]: 
[15]: 



Joshi et alJl2003l): 



McConnachie et al. 



141:lBrownetalJ<2004: 
( 2005i): ri61: lRibasetalJ (2005 ) 



Following the procedure already used in the LMC to ob- 
tain precise distance me asurements to several eclipsing bi- 
nary (EB) systems (see Fitzpa trick et alJ 120031 and refer- 
ences therein), a project was started in 1999 t o obtain a di- 
rect distance determinat ion to M3 1 from EBs jRibas & Jordi 
120031 iRibas et a"fl |2004i). The methodology involves, at least, 
two types of observations: photometry, to obtain the light 
curve, and spectroscopy, to obtain the radial velocity curve 
for each component. By combining the two types of observa- 
tions the individual masses and radii for each of the compo- 
nents can be obtained. The remaining needed parameters (lu- 
minosities, temperatures and line-of-sight extinction) can be 
obtained from the modelling of the spectral energy distribu- 
tion ( see, e.g.. IRibas et alJ l2002) or the modelling of spectral 
lines ( Ri bas et alJl2005l hereafter Paper I). The combination 
of the analysis results yields an accurate determination of the 
distance to the EB system, and hence, to the host galaxy, which 
is the main objective of the project. But in addition to providing 
distance measurements, the resulting stellar physical properties 



can be used as powerful diagnostics for studying the structure 
and evolution of stars that have been born in a chemical envi- 
ronment different from that in the Milky Way. 

The process briefly described above requires high quality 
light curves to obtain precise fundamental properties. By the 
time this project started, the highest quality EB light curves 
were those obtained by the DIRECT group fsee lMacrilEooi 
and references therein), but the scatter was too large for a reli- 
able determination of the elements. Therefore, we began a new 
photometric survey to obtain high quality EB light curves (see 
§0. With the application of the difference image analysis algo- 
rithm in the data reduction process (see ^3), precise photometry 
was obtained not only for the EBs in the field, but also for all 
the variable stars. 

Two photometric catalogs were compiled (see |4). The ref- 
erence catalog contains photometry, in both B and V passbands, 
of 236238 stars from which 3 964 are also in the variable star 
catalog (containing 437 EBs and 416 Cepheids). For some of 
the best quality light curves, corresponding to the brightest EB 
systems - 24 in total, - a preliminary determination of the or- 
bital and physical properties is presented (see §0. These sys- 
tems can be considered potential targets for distance determi- 
nation and, indeed, already one of them wa s used t o derive the 
first direct distance determination to M31 (Paper I). 



2. Observations 

The observations were conducted at the 2.5 m Isaac Newton 
Telescope (INT) in La Palma (Spain). Twenty one nights in five 
observing seasons were granted between 1999 and 2003 (one 
season per year). The Wide Field Camera, with four thinned 
EEV 4128x2148 CCDs, was used as the detector. The field of 
view of the WFC at the INT is 33f8 x 33:8, with 0'.'33 pixel 1 
angular scale. With this configuration, we found that an expo- 
sure time of 15 minutes provides the optimum S/N for stars of 
V ~ 19 mag and ensures that no significant luminosity varia- 
tion occurs during the integration for EBs with a period longer 
than one day. 

All observations were centred on the same field located 
at the North-Eastern part of M31 (a = 00 h 44 m 46 s 5 = 
+41°38'20")- The field of study was selected to overlap with 
the DIRECT fields A, B and C, since the initial main objective 
of the project was to obtain high quality photometry of already 
known EBs. To reduce blending problems in such a crowded 
field (see Fig. all the images with PSF FWHM larger than 
3" were rejected. As a result, 265 and 259 frames, with median 
seeings of 1 '.'3 and 1 '.'2, were selected for further analysis in 
the B and V passbands, respectively. In addition, during two 
nights of the 199 9 observ i ng run and one night of the 2001 
observing run, 1 5 lLandolil dl992h standard stars from 3 fields 
(PG023 1 +05 1, PG1633+099 and MARK A) were observed. In 
Tabled a complete list of all the images used in the data reduc- 
tion process for the five observing seasons is presented. 
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Table 2. Number of images used in the data reduction process 
for each of the five observing seasons. 



Frame 


1999 


2000 


2001 


2002 


2003 


Total 


Granted nights 


3 


4 


4 


5 


5 


21 


Bias 


5 


21 


10 


46 


24 


106 


Darks 











6 





6 


B flat-fields 


20 


31 


19 


21 


42 


133 


V flat-fields 


23 


17 


14 


27 


25 


106 


B Landolt fields 


12 





8 








20 


V Landolt fields 


16 





11 








27 


SM31 field 


41 


45 


45 


54 


80 


265 


VM31 field 


40 


45 


43 


49 


82 


259 


Total 


157 


159 


150 


203 


253 


922 



3. Data reduction 

3.1. Data calibration 

The science images for each of the five observing seasons were 
reduced using the IRAF 1 package. Only bias, dark and fiat-field 
frames corresponding to the same observing season were used 
to reduce a given science image. 

Each WFC image was separated into four overscan- 
subtracted images, one per CCD frame. After this step, each 
CCD frame was treated separately during the whole reduc- 
tion process. Bad pixels 2 were also corrected for all images 
through linear interpolation from the adjacent pixels. Although 
this method is not optimal for reducing photometric data, it 
yields best results for cleaning the calibration images (bias, 
dark and fiat-field images). 

A master bias was built for each observing season and CCD 
frame. Each pixel of the master bias is the median value of the 
same pixel in all the bias frames. The master bias was sub- 
tracted from the 2002 dark images. The dark images of the 
other observing seasons were not used because the dark cur- 
rent was found to be negligible. The bias-subtracted dark im- 
ages were combined in the same way as the bias images (with 
a median) to produce the 2002 master dark frame. 

A master flat-field image was produced for each observing 
season, CCD frame and filter. The bias and, for the 2002 im- 
ages, dark frames were subtracted from each flat-field frame. 
The resulting frames were corrected for non-linearity effects 3 
and averaged to produce the master flat-field image. 

The science images were processed by subtracting the bias 
frame (and also the dark frame for the 2002 images) and cor- 
rected for non-linearity and flat fielding. The resulting images 
are free from most of the instrumental effects, except for a small 
area close to the North-Eastern corner of the field that is asso- 
ciated with field vignetting (see Fig.Q. 

1 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the NSF. 

2 Bad pixel positions were obtained from the webpage: 
http://www.ast.cam.ac.uk/~wfcsur/technical/pipeline/ 

3 Coefficients were obtained from the webpage: 
http://www.ast.cam.ac.uk/~wfcsur/technical/foibles/ 



3.2. Difference image analysis 

Given that the field under study is highly crowde d, a pack- 
age b ased on the image subtraction algorithm dAlard & Luptonl 
1998) was used to obtain the best possible photometric preci- 
sion. This technique has the advantage that variable stars are 
automatically detected and that precise photometry can be ob- 
tained even in highly crowded fields. We have used our own 
implementation of the dif ferenc e image analysis package (DIA) 
developed by IWozniakl |2000). The image subtraction algo- 
rithm requires a high quality reference image. This image was 
created during the first step of the process by the DIA for each 
CCD and filter. The combination of the 15 best seeing images 
produced two reference images, one with 0'.'9 FWHM for V 
(see Fig.Q) and one with 1 '.'0 FWHM for B. 




Fig. 1. Reference image of the selected field of view in V filter. 
The center of the field is ata = 00 h 44 m 46 s and 8 = +41°38'20" 

The main part of the process is the determination of the 
convolution kernels. The kernels are constructed to match the 
reference image with each one of the remaining frames under 
study (supposedly with different seeing values). The DIA ker- 
nels are composed of Gaussian functions with constant widths 
and modified by polynomials. In our case, the kernels consist 
of three Gaussians with polynomials of orders 4, 3 and 2. A 
space-varying polynomial of order 4 has also been used to ac- 
count for the spatial variation across the image. The resulting 
kernels are convolved with the reference image and subtracted 
from the original frames to obtain images containing only ob- 
jects with brightness variations (plus noise). 

The variability image is constructed from the differentiated 
images. This variability image contains all pixels considered to 
have light level variations. A pixel is considered to vary if it 
satisfies at least one of the following two conditions: 
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- There are more than 4 consecutive observations deviating 
more than 3<x (where <x is the standard deviation) from the 
base line flux in the same direction (brighter or fainter). 

- There are more than 10 overall observations deviating more 
than Act from the base line flux in the same direction. 

The variability image is used to identify the variable stars. 
The reference image PSF is correlated with the variability im- 
age and all local maxima with correlation coefficients in excess 
of 0.7 are considered to be variable stars. Once the variable 
stars are located, PSF photometry is performed on all the differ- 
entiated images to obtain the differential fluxes. The differential 
fluxes can be transformed to the usual instrumental magnitudes 
through the following expression: 

nn = C -2.5 logC/b + flA/S) (1) 

where Af is the differential flux of a star in the z'-th observation, 
fo is the base line flux on the reference image, C is a zeropoint 
constant and m, is the instrumental magnitude of a star in the 
;-th observation. 

The photometric package of DIA does not properly account 
for contamination from nearby stars in crowded fields (such as 
our case). Therefore, the fo values obtained with the DIA can 
potentially be severely overestimate d. To obtain p recise values 
for f , DAOPHOT PSF photometry dStetsonll987l) was applied 
to the reference image. Given that DIA and DAOPHOT pho- 
tometry use different PSF definitions and functional forms, a 
scaling factor (a) is neede d to transform the D IA fluxes into the 
DAOPHOT flux units (see IZebrun etal.l200ll for an extensive 
discussion). A synthetic image was created to compute the scal- 
ing factors. The synthetic image contains the representation of 
the DAOPHOT PSFs at the position of each variable star. The 
scaling factor at the position of each variable star was obtained 
from the comparison with the DIA photometry of each PSF. 
For similar PSFs, the scaling factor can be obtained with an 
error of 0.3% or less. 

This process provided instrumental magnitudes for all the 
variable stars as well as all the stars in the reference image. 
Two main sources of potential systematic errors exist with this 
procedure, fo and a, but they affect the amplitudes of the light 
curves. However, as it is shown below in |5] the amplitudes of 
the fitted EBs are perfectly compatible with other results given 
in the literature. 

3.3. Standard photometry 

To tra nsform instrum ental to standard magnitudes, the ob- 
served EandoB 1 19921) standard stars were used. Although the 
standard stars were observed at different times during the night, 
the number of collected frames is not enough to accurately de- 
termine the atmospheric extinction. Therefore, as a first step, 
the M31 images obtained during the same night as the stan- 
dard stars (30 images in total) were used to compute the co- 
efficients needed to account for the atmospheric extinction. 
Aperture photometry of 20 bright, but not saturated, and iso- 
lated stars on the M3 1 frames was performed. For each night 
and filter, a linear fit, with dispersions ranging from 0.008 mag 
to 0.03 mag, was obtained. 



As a second step, aperture photometry was performed on 
the standard stars. The instrumental magnitudes obtained were 
corrected with the atmospheric extinction coefficients. The re- 
sulting values were used to find the transformation coefficients 
between the instrumental magnitudes and the standard magni- 
tudes. For each night and CCD the following relationship was 
determined: 

v-V = A I +A 2 (fi-V) (2) 
b-v = A 3 +A A (B- V) (3) 

where b and v are the instrumental magnitudes, B and V are the 
standard magnitudes and A, are the transformation coefficients. 
The resulting mean standard deviation of the fits is 0.02 mag. 

The transformation coefficients cannot be directly applied 
to the M31 instrumental magnitudes because the latter are 
based on the reference image, which is a combination of differ- 
ent images. In addition, the transformation coefficients are de- 
termined from aperture photometry, while the M3 1 instrumen- 
tal magnitudes have been determined from PSF photometry. 
For this reason, the only way to obtain standard magnitudes for 
the M3 1 stars is to apply the transformation coefficients to the 
frames obtained during the same nights as the standard stars. 

PSF photometry is needed to find precise standard magni- 
tudes of a reasonable number of objects in the M3 1 field since 
aperture photometry can only be applied to isolated stars. A 
sample of bright and isolated stars was employed to determine 
a scaling value used to transform the PSF magnitudes into aper- 
ture photometry values for each one of the 30 M3 1 frames. The 
scaling values and the transformation coefficients were used 
to obtain standard magnitudes for 18 426 objects in the field. 
A systematic difference was observed for the standard magni- 
tudes obtained from the 2001 frames and, therefore, the ob- 
tained values were rejected. The standard magnitudes of every 
star in all the remaining frames (22 in total) were averaged and 
the standard deviation was considered a good estimation of its 
uncertainty. 

The standard magnitudes resulting from this process were 
used to determine new transformation coefficients (A,) for the 
reference images. Only 534 non-variable and non-saturated ob- 
jects detected in all the frames and with an error below 0.04 
mag were used for this purpose. This sample has good color 
and magnitude coverage, with -0.3 < B-V < 1.7 and 
17.5 < V < 21.0, providing fits with dispersions ranging from 
0.013 mag to 0.019 mag, depending on the CCD. The obtained 
coefficients were applied to the reference images. The objects 
detected in the V frame have been cross-matched with the ob- 
jects detected in the B frame. Only the objects identified within 
0'.'33 (one pixel) have been included in the reference catalog, 
providing standard photometry for 236238 objects. The error 
distributions, shown in Fig [2 have 37241 objects with errors 
below 0.1 mag in both passbands. The limiting magnitudes of 
the photometric catalog are around V ^ 25.5 and B 26.0, and 
it is estimated to be complete up to V - 22.3 mag and B 23.5 
mag. The same process was applied to transform the m,- values 
in Eq.[2to standard photometry in both B and V passbands for 
a total of 3 964 variable stars. 
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Fig. 2. Error distribution for the entire photometric catalog con- 
taining 236 238 objects with photometry in both V and B pass- 
bands. 

3.4. Astrometry 

Standard coordinates for all detected objects were computed 
from several reference stars in the GSC 2.2. 1 4 . To ensure that 
the reference stars are uniformly distributed, each CCD was di- 
vided into 3x6 sectors. Three reference stars (with V - 18) 
were identified manually in each sector, yielding an initial 
list of 54 reference stars per CCD. Fits using third order lin- 
ear equations with 2cr scatter clipping were applied to deter- 
mine the transformation of coordinates. After the iterative pro- 
cess at least two reference stars had to remain in each sector. 
Otherwise, an additional reference star was selected in the cor- 
responding sector and the entire process was repeated. 

The resulting coordinates have been compared with those 
in the GSC 2.2. 1 catalog. A total of 724 objects with a position 
difference of less than 3" in the two catalogs were identified. 
As can be seen in Fig. [3] no systematics trends appear in the 
comparisons, which have dispersions of cr a = 0'.'16 and cr 6 = 
0'.'12. 

3.5. Periodicity search 

Since the observations were performed in five annual observing 
seasons, the data is very unevenly sampled. To overcome this 
major drawback for the determination of periodiciti es, a pro- 
gram based on the analysis of variance dSchw arzenber g-Czernvl 
Il996l) was used to compute the periodograms for each variable 
star and filter (B and V). The periodograms were computed 
with 7 harmonics between 1 and 100 days because this setup 
was found to provide optimal results for the detection of EB 
systems, which constitute the main objective of the present sur- 
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Fig. 3. Differences between the coordinates obtained from our 
procedure and the coordinates from the GSC 2.2.1 catalog for 
724 identified stars (circles). The line represents a running av- 
erage. 

vey. The lower limit in the period search is set from the consid- 
eration that EBs with periods below 1 day are likely to be fore- 
ground contact systems and therefore not of the interest of the 
survey. In addition, given the relatively long exposure time for 
each observation in the present survey (15 minutes), important 
brightness variations during the integration time could occur 
for the short period variables. These factors make the compu- 
tational effort of extending the period search below 1 day not 
worthwhile. On the opposite end, the upper limit of 100 days 
was selected to detect as many long period Cepheids as possi- 
ble. 

For each variable star, two periodograms (one in B and 
one in V) were used to obtain a consistent period determi- 
nation. The resulting light curves were visually inspected to 
locate well-defined variability patterns, thus leading to the 
identification of a total of 437 EBs and 416 Cepheids. For sev- 
eral EBs and Cepheids, the period determination was a multiple 
of the true period and, therefore, it had to be recomputed. 

For each variable star, the time of minimum of both light 
curves (B and V) was used to compute a reference time. The 
time of minim um was computed from the Fourier series (F (7) (f) 
in Schwarzenbe rg-Czern\ll996l) that best fit each light curve. 

4. The catalogs 

4.1. Reference catalog 

Photometry of the 236 238 objects detected in the reference im- 
ages was grouped into the reference catalog. The reference cat- 
alog (Table 3) contains the object identifier, the right ascen- 
sion, the declination, the V standard magnitude in the refer- 



6 



F. Vilardell et al.: Eclipsing binaries suitable for distance determination in the Andromeda galaxy 



ence image, the V standard error, the B standard magnitude 
in the reference image, the B standard error and, in case that 
there is a previous identification, the corresponding identifier. 
The catalogs cross-match ed with our referen ce catalog are GSC 
2.2.1 DIRECT, LGGS llMassev et al.ll2006l) and iTodd et ail 
J2005L hereafterlToll). Photographic catalogs are identified in 
the previous works, specially in LGGS, and they are not cross- 
identified here. 

According to the IAU recommendations, the identifier 
was built from the object position and taking into account 
the observational angular resolution. The resulting format is 
M31JHHMMSSss+DDMMSSs. For the variable stars, the 
acronym was changed to M31V to indicate that they are also 
in the variable star catalog (see 34. 2> . 

To show the general photometric properties of the obtained 
catalog, the color-magnitude diagram for the 37 241 objects 
with a photometric error below 0.1 mag in both B and V is 
shown in Fig. 0] From this diagram it can be seen that most 
of the stars in this catalog are stars at the top end of the main 
sequence, with some foreground giants and stars at the tip of 
the red giant branch. 

14 1 1 1 1 — 1 1 1 1 — 1 1 1 1 — 1 1 1 1 — 1 1 1 1 m 1 1 1 1 — m - 

L5- 




-1 -*Ji V U.b I Lb 2 2.3 a 



B-V 

Fig. 4. Color-magnitude diagram for 37 241 objects with pho- 
tometric errors lower than 0. 1 mag in both B and V in the ref- 
erence catalog (grey dots) and 437 EBs (black circles). Stars 
above the dashed lines are saturated. 

The obtained standard photometry depends on a number 
of zeropoint calibrations (§ 33. 21 and 13. 3> . Taking into account 
that an error in the standard magnitudes has a direct impact 
on the distance determination, it is extremely important to en- 
sure that no systematic errors affect the photometry. For this 
reason, it is worthwhile to check the consistency of the stan- 
dard magnitudes given here. Fortunately, there are two CCD 
surveys with B and V photometry of a large nu mber of stars 
that overlap with the field under study: DIRECT dMacril2 004) 



and LGGS dMassev et alJkOOfil) . The non-saturated objects of 
our survey (V > 17.5 and B > 18.0) were cross-identified 
with the magnitudes reported by the DIRECT group 5 and by 
the LGGS group 6 . For the DIRECT survey, a total of 14717 
and 7 499 objects were identified in V and B passbands, re- 
spectively (see Fig. |3- The comparison with the LGGS survey 
provided 36 353 common objects (see Fig.|6jl. 




1.7 La LO 2U 21 22 23 24 



Magnitude 

Fig. 5. Differences with DIRECT standard magnitudes versus 
the magnitudes in the reference catalog. The comparison has 
14717 stars in V and 7 499 in B. The dashed lines show the 
limiting magnitudes. 

Although some trend is observed for the DIRECT B magni- 
tudes, the LGGS B magnitudes are completely consistent with 
the magnitudes in our survey. Therefore, the B standard mag- 
nitudes in the present paper are expected to have no systematic 
error. Regarding the V magnitudes, some low-level systemat- 
ics have been observed with the two catalogs. The DIRECT 
V magnitudes are larger by about 0.02 mag than the values in 
the reference catalog, and the LGGS comparison shows an in- 
creasing trend towards larger magnitudes. The 5 768 objects 
with V magnitudes in the three catalogs were used to study the 
observed systematics. The comparisons between the three cat- 
alogs always reveal some kind of trend and, therefore, it was 
impossible to conclude weather there are any (low-level) sys- 
tematics affecting the V standard magnitudes. In any case, the 
differences for the brightest stars (V < 20.5 mag) are well be- 
low 0.03 mag. Taking into account that all the EBs than can 
be used for distance determination have V < 20.5 ( 35.21 . this is 
the maximum systematic error that can exist in their V standard 
magnitudes. 



5 Data source: ftp://cfa-ftp.harvard.edu/pub/kstanek/DIRECT/ 

6 Data source: http://www.lowell.edu/users/massey/lgsurvey/ 



F. Vilardell et al.: Eclipsing binaries suitable for distance determination in the Andromeda galaxy 



7 



i - 




L7 LS L"J 2U 21 22 23 24 



Magnitude 

Fig. 6. Differences with LGGS standard magnitudes versus the 
magnitudes in the reference catalog for a total amount of 36 353 
stars in each passband. The dashed lines show the limiting mag- 
nitudes. 

4.2. Variable star catalog 

The 3 964 variable stars identified in the reduction process were 
grouped in the variable star catalog (Table 4). For each variable 
star, the corresponding identifier (see 34. 1> . the right ascen- 
sion, the declination, the inte nsity-averaged V magnitude (Eq. 
(9) in lSaha & HoesseJll990l) . the mean V error, the intensity- 
averaged B magnitude, the mean B error, the number of ob- 
servations in V, the number of observations in B, the refer- 
ence time (in HJD), the period (in days) and a label, indicating 
weather the variable star was identified as an EB or a Cepheid, 
are provided. 

Although a period estimate is given for each variable star, 
obviously not all the variable stars are periodic. When com- 
paring the window function with the period distribution for all 
variable stars (see Fig. [7}, it is clear that many of the period de- 
terminations are just an alias introduced by the window func- 
tion, specially for those over five days. Therefore, a significant 
number of the variable stars are, in fact, non-periodic vari- 
ables or variables with a period out of the studied range (1-100 
days). All light curves (time series photometry) are also pro- 
vided (Table 5) together with the variable star catalog. For each 
variable star and observation, the observed time (in JD), the 
standard magnitude (V or B) and the standard error are given. 

5. Eclipsing binaries 

5.1. Complete sample 

The color-magnitude diagram (see Fig. 0} reveals that most 
of the detected EB systems contain high-mass components be- 
longing to the top of the main sequence. Taking into account 



too i | i 




O 0.1 0.3 'J. 3 <JL4 0.5 0.4 0.7 (J.S Ely 1 
Fiajiicjiijf (days I 



Fig. 7. Top: Window function for the studied period range 
(lines) and frequency distribution for all the stars in the variable 
star catalog (histogram). Middle: Frequency distribution for the 
Cepheids in the variable star catalog. Bottom: Frequency dis- 
tribution for the eclipsing binaries in the variable star catalog. 



that most of the systems have periods shorter than 10 days (and 
all of them have periods shorter than 30 days), a large fraction 
of EB systems are non-detached systems. 

Part of the images in the present survey were also reduced 
by|T05 and reported the detection of 127 EB systems. Of those, 
123 stars have been detected in our reference catalog, 92 are in 
the variable star catalog and 90 have also been classified as 
EBs. The maximum V max and B max values for the 90 systems 
identified were matched with the magnitudes presented in lT05l 
(see Fig.|8j. A clear systematic trend can be observed in the V 
filter. In addition, Fig.|8]also reveals that the Bjj^ magnitudes 
are lower than our B max values. 

When inspecting the origin of such discrepancy , we ob- 
served that all the EBs labeled f2BEB and f3BEB in lT05l had 
lower magnitude values than the EBs in our variable star cat- 
alog. Since 14 of these EBs have also DIREC T V max magni- 
tudes , and 5 of them have DIRECT B max values jKaluznv et alJ 
ll998tlStanek et all 199811 19991 iBonanos et alJ2003h . the mean 
differences were computed (Table |6jl. Although the number of 
cross-matched EBs is small, their mean differences seem to 
indicate that the magnitudes reported in T05 for the f2BEB 
and f3BEB EBs suffer from a systematic offset. On the other 
hand, the magnitude values for the remaining EBs (flBEB and 
f4BEB) are fairly compatible in the three catalogs. 

5.2. Light curve fitting 

From all the variable stars identified as EBs only those sys- 
tems with a precise determination of their fundamental proper- 
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18.5 19 19.5 20 20.5 21 21.5 22 22.5 23 
Magnitude at maximum 

Fig. 8. Differences bet ween the DIRECT maximum magni- 
tudes (open squares), the T05 maximum magnitudes (filled cir- 
cles) and the maximum magnitudes for all the crossed EBs of 
the variable star catalog. 

Table 6. Mean differences for the magnitudes at maximum of 



the EBs labeled f2BEB and f3BEB in T05 



ignitui 

Frna 



Comparison 



Differences 



i/mux T/max 

"direct 

|/m;ix ymax 

V T()S DIRECT 



0.062±0.047 
-0.104±0.015 
-0.166±0.057 



nmax nmax 

DIRECT 
£>max _ gmax 

ntTTaT^ nmax 

DIRECT 



0.054±0.039 
-0.214±0.026 
-0.268±0.046 



ties can be used as distance determination targets. Taking into 
account that medium-resolution spectra are needed to obtain 
precise fundamental properties, and considering the fact that 
the largest currently available facilities are the 8-10 m class 
telescopes, only the brightest EB systems (with y mean < 20.5) 
were selected for further study. In addition, the most precise 
fundamental properties can be achieved only for those systems 
with deep eclipses, and thus we further selected only those EBs 
with AB > 0.2 mag and AV > 0.2 mag. 

The criteria above provided a list of 29 systems from the 
initial sample of 437 EBs. Of these, 5 systems were rejected 
because of the large scatter in their light curves (probably 
from the contamination of a brighter nearby star), leaving a 
total amount of 24 EBs selected for detailed further analy- 
sis. The particular properties of each EB can have an impor- 
tant influence on the precise determination of its fundamen- 
tal parameters. To obtain the best distance determination tar- 
gets, a prelimi nary fit was per formed, with the 2003 version of 
the lWilson & Devinnevl(ll97ll) program (W-D), to the selected 
sample of 24 EBs. 



The fitting process was carried out, in an iterative manner, 
by considering both light curves (B and V) simultaneously. For 
those systems previously identified by the DIRECT group their 
V light curve was also included in the fitting process. This pro- 
vided an additional check on the consistency of the resulting 
photometry, because the DIRECT photometry was calculated 
from a PSF rather than using a DIA algorithm. Once a con- 
verging solution was achieved, a Act clipping was performed 
on all the light curves to eliminate (normally a few) outlier ob- 
servations. 

A certain configuration has to be assumed for each solution 
of W-D. Two basic configurations were considered for each 
EB system: detached and semi-detached. Generally, the final 
fit in each one of the configurations provided some clues about 
the real configuration of the EB system. However, all systems 
showing non-zero eccentricity (secondary eclipse phase not ex- 
actly at 0.5) were considered to be detached EB systems. In 
some specific cases, a sinusoidal O - C was observed when 
studying the fit residuals, with one quadrature being brighter 
than the other. This is known as the O'Connell effect. In the 
case of interacting high-mass systems, such effect can be ex- 
plained by the presence of an equatorial hot spot on the sur- 
face of one the components. Such hot spot is supposed to the 
consequence of impacting material arising from mass transfer 
between the components in a Roche-lobe filling configuration 
(semi-detached). The exact parameters of the spot cannot be 
obtained from the ava ilable data (strong degeneracies exist; 
Fitzpatric k et all20 03) and, therefore, the spot parameters pre- 
sented here are not necessarily physically valid but only capa- 
ble of providing a good fit to the variations in the light curves. 

Each EB is a particular case and careful selection of the ad- 
justable parameters was performed individually. In general, the 
parameters fitted with W-D are: the time of minimum (f m ; n ), the 
period (P), the inclination (z), the effective temperature of the 
secondary (Tffl, the normalized potential of primary (Oi), the 
normalized potential of secondary (Qq), for detached systems 
only, and the luminosity of the primary (Lj). For the eccen- 
tric systems, in addition, the eccentricity (e) and the argument 
of the periastron (a>) were also fitted. In some cases of semi- 
detached systems, the mass ratio (q) was treated as a free pa- 
rameter. Finally, a third light {li) contribution was included for 
the final solutions of the EB systems. Taking into account that 
the field under study is extremely crowded, it is not surpris- 
ing that some EBs may have a significant third light contribu- 
tion (see 35.31 for a more extended discussion on this topic). 
In addition, the obtained values of the third light contribution 
provide a first check on the realistic determination of the funda- 
mental properties and give an indication that the scaling factor 
determinations (a in Eq. |2) are correct. Indeed, an error in the 
scaling factor has exactly the same effect as a third light contri- 
bution. The limb-darkening coefficients (square-root law) were 
computed at each iteration from Kurucz ATLAS9 atmosphere 
models and the gravity brightening coefficients, as well as the 
bolometric albedos, were fixed to unity (assuming components 
with radiative envelopes). 

Note that most of the solutions have a large number of free 
parameters. To some extent, the values presented here rely on 
the adoption of the mass ratio and, therefore, the fit results (in 
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Fig. 9. B (up) and V (down) light curves for the 24 fitted EBs. The corresponding Wilson-Devinney fits are also plotted (lines). 
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Table Q and Fig. [5} should be regarded as preliminary until ra- 
dial velocity curves are obtained. For some of the binaries the 
values given in TableQneed additional comments: 

- M31VJ00444528 +4 128000. This is the brightest EB sys- 
tem detected. A converging fit could not be achieved be- 
cause of variations of the light curve at the second quadra- 
ture. The reason is so far unidentified but it could be due 
to a varying hot spot, an accretion disk or intrinsic variabil- 
ity of a supergiant component. This makes the system an 
interesting target for future follow-up observations. 

- M31VJ00442326+4127082. The DIRECT light curve re- 
veals a variation in the argument of periastron thus in- 
dicating the presence of rapid apsidal motion (1.9+0.6 
deg year" 1 ). Even taking into account the important third 
light contribution, the radial velocity curves can provide 
valuable information to study the real nature of this de- 
tached, eccentric system. 

- M31VJ00443799+4129236. Radial velocity curves have 
already been acquired with the 8-m Gemini North telescope 
and mor e accura te fundamental parameters have been de- 
rived, in PaperjJ, from the fit of both radial velocity and 
light curves, as well as from the fit of the atmospheric pa- 
rameters from the spectra. However, we decided to consider 
this system within the current homogeneous analysis and 
thus the comparison of the parameters presented here and 
those in lPaoer 3 provides an idea of the accuracy of the pa- 
rameters in semi-detached systems when the mass ratio is 
unknown. 

- M3 IV J0045 1973+4145048. A likely large error in the 
scaling factor causes a large negative V third light contri- 
bution. In any case, a correct third light contribution is ob- 
tained in the DIRECT light curve and, therefore, the funda- 
mental parameters are expected to be correctly determined. 

- M31VJ00444914+4152526. An equatorial hot spot was 
assumed for this EB system. The best fit yields a spot with 
a temperature ratio of 1.2, at a longitude of 246 + 14 deg 
and a size of 17 + 3 deg. 

- M31VJ00442478+4139022. The large scatter of the 
DIRECT V light curve prevented a reliable fit from being 
achieved. Consequently, the DIRECT light curve was not 
used. 

- M31VJ00445935 +4 130454. A very slight O'Connell ef- 
fect was observed in this EB system, suggesting the pres- 
ence of a hot spot. The best fit, with a temperature ratio of 
1.1, yields a spot at 280 + 60 deg and a size of 2 1 ± 1 1 deg. 

The obtained values provide a good estimation of the po- 
tential of each EB system to become a good distance determi- 
nation target. In general, the best distance determination targets 
should have low third light contribution and a luminosity ratio 
close to unity. On one hand, an important third light contri- 
bution would lead to an underestimation of the observed stan- 
dard magnitudes for the EB system. Although the third light 
can be modelled, the resulting values may have uncertainties 
that are large enough to have a sizeable negative effect on the 
accuracy of the distance determination. On the other hand, EBs 
with luminosity ratios close to unity will have spectra with vis- 
ible lines for the two components, thus making possible the 



determination of accurate fundamental properties. Taking into 
account the obtained values, several EBs have already been se- 
lected as targets to obtain radial velocity measures and, in fact, 
already one of th em has been used to determine the first direct 
distance to M31 (Paper I). 

5.3. Blending 

Special care has to be taken with the magnitudes reported in 
the present paper. As mentioned above, the FWHM of the ref- 
erence i mages is around 1". At the distance of M31 (~770 kpc, 
Paper I) this corresponds to more than 3.5 pc. Therefore, it is 
highly probable that a large fraction of the identified objects 
in the reference catalog are, in fact, a blend of several stars. In 
addition, the blending effect can significantly modify the ob- 
served amplitudes of variable stars an d therefore their derived 
mean magnitudes can also be in error. iMochei ska et alJd2000l) 
studied this problem for Cepheid variables in M3 1 and arrived 
at the conclusion that blending affect the V Cepheid flux at the 
19% level. 

It has been noted before JtoI that such blending would 
seriously hamper the accurate determination of fundamental 
properties and distances of M31 EBs. While this is true for 
Cepheid variables, EB systems are a particular case of vari- 
able stars in which their variations obey a well-defined physical 
model. Therefore, blending effects are easily and seamlessly 
incorporated in the EB modelling as an additional (third) light 
contribution. While special care has to be taken (e.g., check for 
physical consistency between the third light values for differ- 
ent passbands), the presence of a third light does not preclude 
the determination of reliable properties for the components. In 
any case, no systematic errors in the distances of EBs will arise 
from blending because the analysis can detect its effects in a 
natural manner. Note that the mean V third light value for the 
24 fitted systems in our analysis is of 10% and it would be 
around 20% from a rough scaling to single st ars. This is in ex- 
cellen t agreement with the value derived by iMocheiska et alJ 
ll2000h . 

6. Conclusions 

Deep and high-quality time-series photometry has been ob- 
tained for a region in the North-Eastern quadrant of M3 1 . The 
present survey is the deepest photometric survey so far obtained 
in M31. The photometry has been checked for compatibility 
with that from other surveys and, therefore, our resulting cat- 
alog can be used as a reference to extend the time baseline of 
future stellar surveys in M3 1 . 

The study of the variable star population has revealed about 
4000 variables. Further analysis is needed for most of them 
for a proper classification but over 800 variable stars have al- 
ready been identified as EBs and Cepheids. The study of the 
Cepheid variables will be the subject of a forthcoming publi- 
cation JVilardell et alJ2 006). In the current paper we have pre- 
sented the analysis of the EB population, which has resulted 
in an extensive list of EBs suitable for accurate determinations 
of the components' physical properties and their distances. The 
catalog provided here constitutes an excellent masterlist to se- 
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lect systems for further analysis. For example, the resulting 
physical properties will be an important tool to study stellar 
evolution of massive stars in another galaxy and, therefore, in a 
completely independent chemical environment. But, in accor- 
dance with the main goal of the the present survey, the tech- 
nique of determining accurate distances from EBs has been 
demonstrated in Paner I to harbor great potential. Full analysis 
of an additional 5-10 EBs from the sample provided here can 
result in a distance determination to M3 1 with a relative uncer- 
tainty of 2-3% and free from most systematic errors. This will 
represent the most accurate and reliable distance determination 
to this important Local Group galaxy. 
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Table 7. Fundamental properties of the 24 EBs with Wilson-Devinney fits (see text for nomenclature). 



Identifier 


Conf." 


ymax 


P 


i 


e 




R,/a" 


R 2 /a» 


[M31V] 




[mag] 


[days] 


[deg] 




[deg] 






J00444528+4128000 


SD 


18.854 


11.543654+0.000211 


78.0+1.3 


0.0 


- 


0.556+0.009 


0.277+0.002 


J00435515+4124567 


D 


19.164 


6.816191+0.000067 


89.1+4.4 


0.105+0.010 


327.4+ 7.8 


0.160+0.007 


0.281+0.012 


J00442326+4127082 


D 


19.242 


5.752689+0.000037 


87.3+2.4 


0.189+0.019 


51.3+ 4.7 


0.215+0.017 


0.267+0.025 


J0045 1253+4137261 


SD 


19.366 


2.358359+0.000010 


75.4+3.0 


0.0 


- 


0.369+0.023 


0.372+0.007 


J00443799+4129236 


SD 


19.428 


3.549696+0.000012 


83.3+0.6 


0.0 


- 


0.330+0.019 


0.390+0.006 


J00442928+4123013 


SD 


19.475 


3.168969+0.000006 


86.4+2.5 


0.0 


- 


0.329+0.014 


0.403+0.005 


J00450486+4 137291 


SD 


19.536 


3.094681+0.000013 


75.7+3.0 


0.0 


- 


0.316+0.021 


0.390+0.008 


J00451973+4145048 


D 


19.572 


8.130670+0.000087 


87.2+2.4 


0.252+0.027 


109.1+ 1.8 


0.216+0.014 


0.149+0.012 


10045005 1+413 1393 


D 


19.756 


5.211976+0.000044 


81.7+1.2 


0.040+0.015 


120.0+14.1 


0.128+0.005 


0.349+0.011 


J004449 14+4 152526 


SD 


19.772 


2.626992+0.000010 


72.5+1.5 


0.0 




0.340+0.015 


0.384+0.005 


J00442478+4139022 


D 


20.034 


4.762310+0.000068 


79.9+1.2 


0.027+0.008 


188.8+33.4 


0.105+0.003 


0.333+0.014 


J00443610+4129194 


SD 


20.038 


2.048644+0.000006 


70.6+1.7 


0.0 




0.290+0.018 


0.369+0.007 


J00425907 +4136527 


D 


20.060 


5.874724+0.000069 


81.0+1.1 


0.163+0.016 


70.8+ 2.2 


0.344+0.013 


0.157+0.011 


J00452554+4145035 


SD 


20.090 


5.009412+0.000047 


82.3+3.4 


0.0 




0.247+0.020 


0.393+0.008 


J00442722+4 136082 


SD 


20.098 


4.518795+0.000016 


86.1 + 1.7 


0.0 




0.308+0.018 


0.344+0.009 


J00453244+4 147425 


SD 


20.141 


2.787856+0.000010 


78.4+1.7 


0.0 




0.238+0.012 


0.342+0.007 


J00445935+4130454 


SD 


20.166 


2.668419+0.000023 


68.1 + 1.9 


0.0 




0.277+0.025 


0.347+0.013 


J00450537+4 133402 


SD 


20.326 


1.769903+0.000005 


75.4+2.5 


0.0 




0.351+0.025 


0.350+0.010 


J00432610+4142010 


SD 


20.340 


4.274429+0.000065 


75.9+3.6 


0.0 




0.206+0.027 


0.413+0.014 


J00445682+4131109 


D 


20.360 


4.207679+0.000036 


87.2+4.5 


0.0 




0.345+0.015 


0.158+0.015 


J00435082+4121517 


D 


20.367 


2.176672+0.000011 


79.6+4.1 


0.0 




0.305+0.026 


0.358+0.031 


J00425743+4137114 


SD 


20.388 


1.916302+0.000008 


73.9+2.9 


0.0 




0.364+0.019 


0.379+0.006 


J00425929+4138169 


D 


20.450 


5.591515+0.000045 


85.5+4.0 


0.0 




0.340+0.011 


0.264+0.016 


J0044501 1+4128069 


D 


20.456 


2.861046+0.000009 


80.2+2.4 


0.0 




0.347+0.013 


0.286+0.013 



"Configuration. (D): Detached, (SD): Semi-detached. 

b R l /a and R 2 /a are the relative radii of the components, where a is the semi-major axis of the system. 



Identifier 
[M31V] 


Tf/Tf 


(Ll/Li)" 


(Li/Li)" 




'3 


•I 


h 


J00444528+4128000 


0.453+0.033 


0.0513+0.0003 


0.0575+0.0004 


0.0575+0.0004 


0.0 


0.0 


0.0 


J00435515+4124567 


0.911+0.018 


2.69+0.09 


2.70+0.08 




0.055+0.032 


-0.005+0.030 




J00442326+4127082 


0.890+0.027 


1.25+0.13 


1.27+0.14 


1.27+0.15 


0.329+0.039 


0.336+0.043 


0.328+0.051 


J00451253+4137261 


0.885+0.020 


0.80+0.08 


0.82+0.09 


0.82+0.12 


0.285+0.056 


0.308+0.059 


0.325+0.070 


J00443799+4 129236 


0.772+0.012 


0.88+0.04 


0.90+0.04 


0.90+0.04 


0.0 


0.0 


0.0 


J00442928+4123013 


0.857+0.011 


1.17+0.05 


1.19+0.05 


1.19+0.06 


0.042+0.035 


0.027+0.032 


0.022+0.034 


J00450486+4137291 


0.787+0.015 


1.01+0.12 


1.03+0.12 


1.03+0.13 


0.016+0.098 


-0.041+0.093 


-0.033+0.096 


J00451973+4145048 


0.894+0.016 


0.34+0.06 


0.35+0.06 


0.34+0.07 


0.041+0.077 


-0.136+0.096 


0.115+0.096 


J00450051+4131393 


0.993+0.005 


7.67+0.35 


7.74+0.38 


7.73+1.47 


0.063+0.025 


0.069+0.024 


0.055+0.086 


J00444914+4152526 


0.957+0.026 


1.20+0.09 


1.21+0.09 




0.003+0.046 


0.009+0.047 




J00442478+4139022 


0.848+0.004 


7.48+0.34 


7.65+0.35 




0.0 


0.0 




J00443610+4129194 


0.946+0.020 


1.51+0.18 


1.52+0.18 


1.52+0.22 


0.067+0.034 


0.051+0.032 


0.016+0.047 


J00425907+4 136527 


1.046+0.021 


0.22+0.02 


0.22+0.02 


0.22+0.02 


0.053+0.017 


0.051+0.029 


0.034+0.023 


J00452554+4145035 


0.810+0.018 


1.72+0.13 


1.78+0.13 


1.78+0.27 


0.169+0.014 


0.162+0.023 


0.227+0.029 


J00442722+4 136082 


0.747+0.008 


0.72+0.05 


0.76+0.05 


0.76+0.06 


0.109+0.009 


0.130+0.010 


0.119+0.014 


J00453244+4147425 


0.728+0.009 


1.14+0.11 


1.20+0.12 


1.20+0.17 


0.036+0.022 


0.018+0.026 


0.093+0.090 


J00445935+4130454 


0.857+0.063 


1.22+0.24 


1.25+0.25 




0.169+0.079 


0.219+0.082 




J00450537+4 133402 


0.940+0.015 


0.88+0.10 


0.89+0.10 


0.89+0.13 


-0.055+0.110 


-0.027+0.102 


-0.037+0.117 


J00432610+4142010 


0.838+0.028 


3.15+0.77 


3.23+0.79 




0.243+0.043 


0.309+0.061 




J00445682+4131109 


0.891+0.024 


0.18+0.03 


0.18+0.03 


0.18+0.04 


0.118+0.076 


0.129+0.087 


0.153+0.117 


J00435082+4121517 


0.789+0.019 


0.89+0.18 


0.93+0.19 




0.245+0.074 


0.211+0.091 




J00425743+4137114 


0.921+0.014 


0.94+0.11 


0.96+0.11 




0.165+0.063 


0.219+0.070 




J00425929+4138169 


0.548+0.011 


0.17+0.02 


0.19+0.02 


0.19+0.02 


0.195+0.042 


0.203+0.049 


0.139+0.045 


J0044501 1+4128069 


0.603+0.016 


0.25+0.03 


0.26+0.03 


0.26+0.04 


-0.012+0.090 


-0.004+0.080 


0.001+0.105 



